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Abstract 

We study the dynamics of N point vortices on a rotating sphere. The Hamilto- 
nian system becomes infinite dimensional due to the non-uniform background 
vorticity coming from the Coriolis force. We prove that a relative equilibrium 
formed of latitudinal rings of identical vortices for the non-rotating sphere per- 
sists to be a relative equilibrium when the sphere rotates. We study the nonlin- 
ear stability of a polygon of planar point vortices on a rotating plane in order to 
approximate the corresponding relative equilibrium on the rotating sphere when 
the ring is close to the pole. We then perform the same study for geostrophic 
vortices. To end, we compare our results to the observations on the southern 
hemisphere atmospheric circulation. 
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1 Introduction 

The interest of studying point vortices on a rotating sphere is clearly geophys- 
ical. This may permit also to understand the motion of concentrated regions 
of vorticity on the surface of planets such as Jupiter |DL93| . The literature is 
now numerous on point vortices on a non-rotating sphere |B771 IKN981 [KN99 
IkJMOOI INOOl IPM98I iLMkOll IBCOII ITEM] b ut only few papers consider a ro- 
tating sphere [F75l 15771 lB85l IKR891 IDP98j . In |F75], the interaction of three 
identical point vortices equally spaced on the same latitude is investigated via 
the /3-plane approximation. It is shown — under some additional assumptions 

— that this configuration is an equilibrium, and that its linear stability depends 
on the strength of the vortices: linearly stable for a negative or a strongly posi- 
tive strength, linearly unstable otherwise. In |B77j . the equations of motion on 
a rotating sphere are given, while in B85| the motion of a single point vortex 
is given. It appears that a single vortex moves westward and northward as a 
hurricane does. In KRJ39], the approach is completely different from ours and 

B77 ; they proved the existence of relative equilibria formed of two or three vor- 
tices of possibly non- identical vorticities. In |DP98j . they model the background 
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vorticity coming from the rotation by latitudinal strips of constant vorticitics 
and they study the motion of a vortex pair (point-vortex pair as well as patch- 
vortex pair). A vortex pair is a solution formed of two vortices with opposite 
sign vorticities rotating around the North pole. A vortex pair moving eastward 
or strongly westward is stable, unless the vortex sizes are too large. Different 
types of instabilities are described for weak westward pairs. 

In this paper, we first recall some basics notions of geometric mechanics in 
Section [3 In particular, the notion of relative equilibrium is defined. In the 
cases of that paper, a relative equilibrium corresponds to a rigid rotation of N 
point vortices about some axis. We then give the equations of motion for point 
vortices on a rotating sphere |B77| in Section [3] The Hamiltonian system is 
infinite dimensional due the background vorticity coming from the rotation of 
the sphere. 

We prove that a relative equilibrium formed of latitudinal rings of the non- 
rotating system persists when the sphere rotates. From LMR01 and |LP02j . we 
know that the following arrangements of latitudinal rings are relative equilibria 
of the non-rotating system: CNv{R,k p p), Djv/i(2-R, k p p), D^d{R, R' ' ,k p p), and 
D2Nh(R e ,k p p) (k p is the number of polar vortices). See Figures^!!! and|3] 




(a) C 4v (R) (b) C 4v (R,p) 



Figure 1: The C nv (R, k v p) relative equilibria with k p = 0, 1. 

In Section^] we give the stability of relative equilibria Cnv(R) and Cn v (R,p) 
for three different approximations or limiting cases: point vortices on a rotating 
plane, geostrophic vortices, and point vortices on a non-rotating sphere. Indeed 
in that particular cases the system becomes finite dimensional and we can use 
the techniques of Section[2]to obtain both nonlinear and linear stability results. 
The stability is determined with a block diagonalization version of the Energy- 
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D ih (2R) D 4d (R,R') 
Figure 2: Relative equilibria D 4 h(2R) and D 4 d{R,R')- 




Figure 3: Equatorial (±)ring D2Nh(R e )- 
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momentum method jU"R99 LP02, and the Lyapunov stability results are modulo 
SO(2). In particular, we compute the (nonlinear) stability of a polygon formed 
of N identical point vortices in the plane together with a central vortex of 
arbitrary vorticity in Appendix B. Our results differ from those of CS99] but 
agree with the linear study of |MS71| . We also improve some stability results 
on geostrophic vortices of MS71| proving that some linear stable configurations 
are actually Lyapunov stable. 

The paper ends with a discussion on the Southern Hemisphere Circulation 
and its relationship with vortices. 

2 Geometric mechanics 

We recall thereafter some basics notions of geometric mechanics. We refer to 
| MR941 IOr98l for further details. 

Let G be a connected Lie group acting smoothly on a symplectic manifold 
(P, w). Consider an Hamiltonian dynamical system (V, vj, H) with momentum 
map J : V — » q* such that the Hamiltonian vector field Xh and the momentum 
map are G-equivariant. A point x e € V is called a relative equilibrium if for all 
t there exists gt G G such that x e (t) ~ gt • x e , where x e (t) is the dynamic orbit 
of Xh with x e (0) — x e . In other words, the trajectory is contained in a single 
group orbit. A relative equilibrium x e £ V is a critical point of the augmented 
Hamiltonian: 

Hz(x) = H(x)-(J(x)-tM,Z) 

for some ( £ g. The vector £ is unique if the action of G is locally free, and is 
called the angular velocity of x e . 

In the case of point vortices on a rotating sphere, we will have G = SO(2). 
Hence relative equilibria are rigid rotations with angular velocity £ € so (2) ~ K 
around the axis of rotation of the sphere. 

The Hamiltonian is G-invariant, but may have additional symmetries. De- 
note by G the group of symmetries of the Hamiltonian. For example, in the case 
of N identical point vortices on a non-rotating sphere, one has G = 5*0(3) x Sjy 
and G = 0(3) x Sjy LMR01 . Let the fixed point set of a subgroup K of G be: 

Fix{K,V) = {x e V | g ■ x = x,Wg £ K}. 

The following theorem permits to determine relative equilibria, it is a corollary 
of the Principle of Symmetric Criticality of Palais |P79| . 

Theorem 2.1 Let K be a subgroup of G, x S Fix(if, V) and /i = J{x). Assume 
that G is compact. If x is an isolated point in Fix(ii', V) PI J _1 (/i), then x is a 
relative equilibrium. 

Note that this result depends only on the symmetries, the phase space and the 
momentum map, and not depends on the form of the Hamiltonian. A relative 
equilibrium obtained with that result is said to be a large symmetry relative 
equilibrium since its isotropy subgroup must be large. 
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To compute the stability of relative equilibria, we use the Energy momen- 
tum method: let x e be a relative equilibrium, fi = J(x e ) € Q* and £ be its 
angular velocity. The method consists first to determine the symplectic slice 

M = T Xc {G„ ■ x,) 1 - n Ker DJ(x e ) 

where 

G M = {g e G | Coad g ^ = /i}. 

The second step consists to examine the definitness of d 2 H^\j^(x e ), and apply 
the following result |Pa92l IUR99| which holds in particular if G is compact: 
If d 2 H^\j^f{x e ) is definite, then x e is Lyapunov stable modulo G^. 

In Section 14.11 we will consider vortices in the plane (point vortices and 
geostrophic vortices), the symmetry group G is SE{2) or SO (2) depending on 
whether the plane is rotating. However we will forget translational symmetries 
since SE(2) is not compact, hence we take G — SO (2) and the previous result 
holds. Moreover we have G M = 50(2) for all fi e so (2)* ~ E. In Section 
14.21 the sphere is non-rotating, thus the symmetry group is G = SO (3) LP02 
LMRj. We have 50(3) M = SO (2) for fx + 0, and 50(3)^=0 = 50(3). Since 
we will consider only relative equilibria with a non-zero momentum ji for that 
section, Lyapunov stable will mean Lyapunov stable modulo SO (2) throughout 
that paper. 

The linear stability is investigated calculating the eigenvalues of the lin- 
earization in the symplectic slice, that is of L_\f = fi^ 1 d 2 H^\^j-{x e ) where VtKj- 
is the matrix of 

The symplectic slice N is a G Xe -invariant subspace. Hence we can perform an 
G x<! -isotypic decomposition of A/", this permits to block diagonalize the matrices 
d 2 H^\f^(x e ) and La/, their eigenvalues are then easier to compute and we can 
conclude about both Lyapunov and linear stability. A basis of the symplectic 
slice in which these matrices block diagonalize is called a symmetry adapted 
basis. The symmetry adapted bases do not depend on the particular form of 
the system, they depend only on the symmetries of the system. 

The different steps of the method are widely detailed in |LP02| which is a 
study on point vortices on a non-rotating sphere. 

A relative equilibrium x e is said to be elliptic if it is spectrally stable with 
d H^\jsf(x e ) not definite. An elliptic relative equilibrium may be Lyapunov 
stable, but this can not be proved via the Energy-momentum method (but 
KAM theory may work). Note also that an elliptic relative equilibrium becomes 
linearly unstable when some dissipation is added to the system DR02 . 

3 Equations of motion and relative equilibria 

In this section we consider N point vortices on a unit sphere rotating with a 
constant angular velocity Q around the axis (Oz). Hence the 50(3) symmetry 
of the non-rotating system breaks to 50(2). The Coriolis force induces a con- 
tinuous vorticity ujq on the sphere: u = luq + loq where loq — \iS(x — Xi), 
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Xi is the vorticity of the vortex Xi, and cun(t = 0) = 2£lcos#. The continuous 
vorticity is not uniform and thus interacts a priori with the singular vorticity 
i.e. the vortices. This interaction makes u>n a function of time. For example, 
without any vortices ui = ojq = 2£lcos# is a steady solution. 

The stream function ifj satisfies Aip = u>. Hence ip — G * u> where G is the 
Green function on the sphere G(x.x') = l/(47r) ln(l — x ■ x'), and we obtain the 
following expression for the stream function: 



1 N 

ip = — } Xi ln(l - x ■ Xi) + G * ton. 

47T * — ' 



The equations governing the motion of the vortices are therefore [B 

N 



E sin^sin(^- - fc) 1 dG * luq 
1 — Xi ■ Xj "-• <™>-- 



N 

, sm Ui cos tt~ — 

A . 

1 — Xi ■ Xi d9i 



. ■ x sin Bicosdj —sin 9 jCDsdiCoa(<f>i — fa) dG * lu^ 
S m9i<f>i = - ^ Xj i- Xi . x . IT 



for alii = 1, . . . , N. Moreover, the vorticity satisfies dtU> + (u-V)u> = (Euler), u 
being the velocity, and this equation can be written equivalently with a Poisson 
bracket: 

w = W,ip} 

where 

df dg df dg 

if, 9} 



d cos 9 dp dtp d cos 9 

for two smooth functions /, g on the sphere. The full dynamical system is 
therefore: 

• ^ sinfljsin^-^) 1 dG * Uq 

t>i = > Aj h -— — [pi,(pi) 

. ^r!, . 1-Xi-Xj sm 9 { dpi 

. „ ; . sin 9i cos 9j - sin 9 a cos 9i cost fa - 0,) <9G * Ufi. . 
sm9i<t>i = - ^ A, 



1 - x l • x,- d9 l 

3 = ^,3=Fi 

tin = {un,ip} 
since the strengths of the vortices are constant. 

Remark. The vorticity must satisfy J uj dS — from Stoke's theorem. It 
should be noted that when the sum A = ^ is non-zero, then the stream 
function remains as before taking ujq = — A/(47r) + J^. Xi5(x — Xi). Since loq 
satisfies J uiq dS — 0, it follows that J oj dS = when wq = 20 cos 9. 
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It is easy to verify that the following quantity (total kinetic energy) is con- 
served and serves as a Hamiltonian: 



Let SO(2) be the group of rotations with as axis the axis of rotation of the 
sphere, and consider the diagonal action of SO(2) on the product of N spheres. 
Clearly the continuous vorticity satisfies 

un(g ■ 6, g ■ 4>, g ■ &i, g ■ <f>i, . . . , g ■ 9n , 9 ■ <I>n) = un(9,<i>,8i,<t>i, ■ ■ -,9 N ,cj) N ) 

for all g e 5*0(2). It follows that H is S'0(2)-invariant and the dynamical 
system is S0{2)-eqm variant. 

Due to the continuous vorticity, the phase space V becomes infinite dimen- 
sional: 



where F{S 2 ) is the set of smooth functions on the sphere. 

The following vector (momentum vector) is conserved and thus provides 
three conserved quantities |B77| : 



The term ^ Aj-Xj corresponds to the moment map for point vortices on a non- 
rotating sphere LMROI]. Indeed in the case of the non-rotating sphere, the 
dynamical system is 5 l O(3)-equivariant, and this leads to three conserved quan- 
tities since SO (3) is of dimension three. In the case of the rotating sphere the 
system is only SO(2)-equivariant, hence the symmetries provide only one con- 
served quantity (which is the z-componcnt of M) . Actually, the three conserved 
quantities come from a general property of 2D incompressible and inviscid fluid 
flows on compact and simply connected surfaces which states that the vector 
J ui-x dS is conserved, no matter the symmetries of the surface we have [B]. That 
vector and the momentum map simply coincide in the case of a non-rotating 
sphere. 

We would like to know if the relative equilibria found on the non-rotating 
sphere persist when the sphere rotates. We call an N-ring a latitudinal regular 



H 





V = {(x u ...,x N ) e S 2 x ••• x S 2 | Xi ±xj if i ^ j} x T{S 2 ) 
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polygon formed of N identical vortices. The following theorem show that rel- 
ative equilibria formed of ./V-rings (together with possibly some polar vortices) 
persist. 

Theorem 3.1 Let x e be a relative equilibrium of the non-rotating system, with 
angular velocity £o> formed of m N m -rings (N m > 2) together with possibly 
some polar vortices. Then (x e ,con(t) = 2f2cos(9) is a relative equilibrium of the 
rotating system with angular velocity £ = £q + fl. 

Corollary 3.2 The relative equilibria 

C Nv (R,k p p), D Nh {2R,k p p), D Nd (R,R',k p p), D 2 Nh(Re,k p p) 

persist when the sphere rotates. 



Proof. Let (a; e ,£o) be a relative equilibrium of the non-rotating system 
formed of k r iV^-rings, a polar vortex being considered as a 1-ring. Hence 

• 1 dG * uon . 



smt 



i) 



1 dG * wjj 







sin $i dOi 

for all i = 1, . . . , N. 

Let the continuous vorticity be 0Jn(t) = 2f2cos# for all time t. We have 

)d(f> does not depend on tfii, and 



= since jG(9,,cj ) 



dG * uj n 
d9 r 



47r J 6=0 



sin 29 



cos 9 sin 9j — sin 9 cos 0; cos( 



i) 



=o 1 — cos 9 cos 9i — sin 9 sin 



n 

2 Je=o 



sin 29 



cos u sin f 



1 - 



9i cos(0- fa) 
1 — cos 9 cos 9i 



d<t>d9 
d9. 



| cos a — cos t)i | sin f 

And we obtain after some calculus that dG g^ n = —rising. 
Thus we proved that 9\ = and ^ = £o + ^ for all z = 1 



cos f 



cosf 



then to prove that ton — 0. One has ujq — 
it follows that 



2m 



. . . , N. It remains 
Since = 0, 



dip kr Nk 
Y fe=i i=i 



sin sin 9k sin(0 — 



cos 9 cos — sin 9 sin 6^ cos(</> — tpj.k) 



where and 9k are respectively the vorticity and the co-latitude of the ring 
fc, and 4>j t k = £fe + 2ir(j — l)/N. It can be shown that this sum vanishes (see 
Appendix A), hence uin — and (x e , o>n(i) = 2f2cos#) is a relative equilibrium 
of the rotating system with angular velocity £ — £o + *3 



Remark. In order to prove existence of relative equilibria, we may think to 
use the Principle of Symmetric Criticality of Section |2 The Principle holds 
since SO(2) is compact, however the fixed point sets Fix(CN,V) are infinite 
dimensional, it is therefore quite a task to find relative equilibria with that 
method. 
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4 Stability of relative equilibria on the rotating 
sphere 

In this section, we unfortunately do not compute the stability of the relative 
equilibria determined in the previous section. Indeed the Hamiltonian system 
is here infinite dimensional and the method described in Section [21 work only 
for finite dimensional Hamiltonian systems. Hence we will give in this section 
the stability of different approximations or limiting cases of the "rotating sphere 
problem" which lead to a finite dimensional Hamiltonian system. 

4.1 Vortices on a rotating plane 

We consider point vortices on a rotating plane in order to approximate the 
Cnv(R) and Cnv{R,p) relative equilibria on a rotating sphere when the ring 
is close to the pole. This approximation can be done in two different manners: 
the "classical" point vortices on a rotating plane, and the "geostrophic" vor- 
tices |S43| . The equivalent of the Cnv (R, p) arrangement on the plane is the 
arrangement C'n (R, p) : a regular polygon formed of N vortices of strength 1 to- 
gether with a central vortex of strength A. When there is no central vortex, the 
configuration is denoted Cn(R) and corresponds to the Cnv{R) arrangement 
on the sphere (see Figure [Q. We assume that N > 3 since the cases N = 2 are 
degenerate due the colinearity of the arrangements. 




(a) d(R) (b) C 4 (R,p) 



Figure 4: The C/v(i?, fc p p) planar relative equilibria. 



Point vortices on a rotating plane. 

Consider a system of point vortices on a plane rotating with a constant angular 
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velocity SI around its normal axis. Contrary to the rotating sphere, the con- 
tinuous vorticity induced by the rotation is here uniform, thus the continuous 
vorticity does not interact with the point vortices. The dynamics of point vor- 
tices on the rotating plane is therefore similar to that for the non-rotating plane 
as the following proposition shows. 

Denote a relative equilibrium x e = (zi,...,Zn) € C N with angular velocity 
£ by (a; e iO an d se t the origin of the plane to be the centre of the rotation. 
Recall that A is the sum of the vorticities, that is A = J^. Aj. 

Theorem 4.1 If(x e ,£) is a relative equilibrium of the non-rotating plane such 
that J2^j z j — an d A ^ 0, then (x e ,£ + Si) is a relative equilibrium of the 
rotating plane. Moreover (x e ,£ + SI) satisfies the stability properties of (x e ,£) 
(as elements ofV/SO(2)). 

Proof. The continuous vorticity induced by the rotation of the plane is 
u>n(x,y) — 2SI, thus the distribution of vorticity for N point vortices Zj = 
pj exp(i<f>j), j = 1, . . . , N is given by u = 2SI + J2j A? ~~ z j)- ^ follows that 
the Hamiltonian is 

where Hq is the Hamiltonian for the non-rotating plane (see Appendix B). 
Since Xjpj = and Xj<j)j = j-§^-, the dynamical system is given by 

{i)-{iy(-m\ 

Since XjZj = and A^O, one can show that the relative equilibrium of the 
non- rotating plane (x e ,£) satisfies 

<^£ = 1 dH ° _ g 
'><:>, ' Xjpj dpj 

for all j (A 7^ is required to insure the existence of the "barycentre" ^ Xj Zj/K). 
Hence (x e , £) with x e = x e and £ = Sl + ^ is a relative equilibrium of the rotating 
plane. 

If a vortex Zj of the relative equilibrium x e is at the origin (a central vortex), 
this discussion is not valid due to the degeneracy of polar coordinates at the 
origin, we then use cartesian coordinates and find the same conclusions. 

We are interested in stability in V/SO{2) (stability modulo SO (2)) where 
SO (2) is the group of rotations around the origin. The "rotating" dynamical 
system is S'0(2)-cquivariant. The momentum map coming from the SO (2) 
symmetry is J = J2j Xjp 2 /2 €R~ so (2)* (see Appendix B), thus H = H a + SlJ. 
Hence 

d 2 H^ e ) = d 2 H(x e ) - {Sl + i)d 2 J{x e ) = d 2 H Q (x e ) - id 2 J{x e ) = d 2 H 0i (x e ). 
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Moreover, the symplectic slice is the same for both the rotating and the non- 
rotating cases, since the stability is investigated in V / SO{2) (and not in V/SE(2) 
as we could do for the non-rotating case). The stability statement follows. □ 

Remark. The previous result is actually due to the fact that the Hamiltonian 
H is perturbed by a function J n (x) — (J(x),Cl) where J is a momentum map 
of a symmetry group G of Hq, and 17 6 q. The Hamiltonian J fi is a collective 
Hamiltonian (that is a function of J) , see |GS84| for a general study of collective 
Hamiltonians. 

Consequently, all relative equilibria of the non-rotating plane formed of reg- 
ular polygons persist to be relative equilibria when the plane rotates, provided 
that the polygons are cocentric with centre of the rotation. In particular, the 
Cn(R) and Cjsr(R,p) relative equilibria persist (the centre of the ring being 
the centre of the rotation). For relative equilibria formed of several regular 
polygons, we refer to |LR96| in which a stability study is also given for some 
particular configurations. 

We can now give stability results for the Cn (R) and Cm (R, p) relative equi- 
libria thanks to the stability study of CV(i?) relative equilibria |CS99| on the 
non-rotating plane and a theorem proved in Appendix B. 

Proposition 4.2 The CV(i?) relative equilibria on a rotating plane are Lya- 
punov stable if N < 7 and linearly unstable otherwise. 

Let En = +1 if N is even, and En = if N is odd. 

Theorem 4.3 The Cn{R,p) relative equilibria with N > 4 on a plane (rotating 
or not) are Lyapunov stable «/max(0, (N 2 -8N + 7 + e N ) /16) < A < (N-l) 2 /4, 
elliptic if (N 2 — 8N + 7 + 6n)/16 < A < 0, and linearly unstable if the previous 
conditions are not satisfied. The C3 (R, p) relative equilibria are Lyapunov stable 
ifO < A < 1, elliptic if X < 0, and linearly unstable if X > 1. 

We have (N 2 — 8N + 7 + e N )/16 < if and only if N < 6, hence there 
exists some linearly stable configurations with a negative central vorticity only 
for N < 6. 



Proof. The stability results of Cn(R,p) relative equilibria on a non-rotating 
plane are given in Appendix B since the proof is quite long. These results will 
hold for a rotating plane by Theorem 14. II □ 

Remark. The Lyapunov stability of Cf(R) (the Thomson heptagon) is not 
proved with an energetic method. Indeed the Hessian is not definite at this 
relative equilibrium, one needs to go to the fourth order to determine the Lya- 
punov stability |M78llCS99| . 
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Geostrophic vortices. 

In his paper |S43| , Stewart considered a single layer of fluid of constant density 
on a rotating disc in order to model the atmosphere of the Earth. Let f2 be the 
constant angular velocity of the disc, z be the complex coordinate on the rotating 
disc, ho be the uniform depth of the fluid at rest and g be the acceleration due to 
gravity. By means of the geostrophic wind equations, he found that the stream 
function corresponding to a vortex at zq with vorticity A is 

* = XK (k\z- z \) 

where n — 2fl/^/gho and Kq is the Bessel function of the second kind. Such a 
vortex is called a geostrophic (or Bessel) vortex. The Hamiltonian for a system 
of N geostrophic vortices z\, . . . , zn with vorticities Ai, . . . , Ajv is 

H =} \i\jKo(n\zi - zj\). 

i<j 

Note that a geostrophic vortex is almost a planar k~ -Euler vortex, the param- 
eter k _1 is the horizontal length scale of the geostrophic vortex. In the case of 
the Earth k _1 is about 2000 kms. If we rescale the radius of the Earth to be 1 as 
the spheres previously considered, K is about 3. Note also that the limiting case 
K = corresponds to logarithmic vortices, that is to classical point vortices. 

Since the Hamiltonian for the geostrophic vortices has the same symmetries 
as the Hamiltonian for the planar point vortices, we follow Section|21to determine 
large symmetry relative equilibria and study their stability: it is easy to show 
that CW (R) and C n (i?, p) configurations are large symmetry relative equilibria 
for the system of planar point vortices, and to find the symmetry adapted basis. 
The symmetry adapted basis for Cn(R,p) is given in Appendix B, while the 
one for Cn(R) is easily deduced from the former. These bases are similar to 
those for the Cnv{R) and Cnv{R,p) configurations of point vortices on a sphere 
|LMR| . 

We set the value of the vorticity of the ring to be 1. The stability results are 
given in Proposition ^. 4l for Cn(R) geostrophic relative equilibria and in Figure 
[5] for Cn(R,p) geostrophic relative equilibria. The linear stability results were 
already found by Morikawa and Swenson ,\ I S 7 1 1 by means of a numerical com- 
putation of the eigenvalues of the linearization. Here we compute numerically 
the eigenvalues of the diagonal blocks of the Hessian 1 and the linearization, and 
thus conclude about both formal (hence Lyapunov) and linear stability. The 
results of this paper agree with their results except for the right stability fron- 
tiers where they differ slightly. This is probably due to a lack of accuracy. The 
advantage of the method described here is that we know that both the Hessian 
and the linearization block diagonalizc, hence it is not necessary to compute the 
nil blocks (gain of time), and we know that the components in it are all zero and 
this contributes to improve the accuracy of the eigenvalues (indeed, if we com- 
pute numerically the components of the nil blocks, we find for example 10~ 13 

1 But the Hessian is obtained analytically. 
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instead of zero, and this leads finally to a lack of accuracy in the computations 
of the eigenvalues). 

Proposition 4.4 The stability of the geostrophic Cjv(-R) relative equilibria de- 
pends on N and k as follows: 

• N > 7: linearly unstable for all n > 0; 

• N = 6: Lyapunov stable for k G [0, 1.28], linearly unstable otherwise; 

• N — 5: Lyapunov stable for k 6 [0,3.75], elliptic otherwise; 

• N — 3, 4: Lyapunov stable for all k > 0. 

Instability on the left part of the figures is due to the ring, while instability in 
right part is due to the central vortex. For N — 3, the linear stability frontier of 
the central vortex (the right one) is almost a vertical line for k € [0, 5], that is K 
in that range has almost no influence on the linear stability of the central vortex. 
For TV = 4, . . . , 7, the stability frontier of the ring (the left one) has the same 
property. On the other hand, the stability frontiers of the central vortex (the 
two on the right) for N > 4 depend on k. Note that the notion of "dependence" 
is linked to the scale of the figure for N = 5, that is (k, A) S [0, 5] x [—2, 10]. 

For N < 7, almost every relative equilibrium stable for k = is stable for 
k G [0, 5] . Thus if a point vortices Cn{R,p) relative equilibrium is stable, then 
so does the corresponding geostrophic relative equilibria. The point vortices 
system is therefore a good model to determine the stability (not instability) for 
N < 7. This becomes less and less true as iV increases. Indeed the greater is 
TV, the more the stability depends on k. 

4.2 Stability results for Q = 

In this section, we give the stability results for the Cnv(R) and Cjy v (R,p) 
relative equilibria on the non-rotating sphere. 

The following theorem was proved first in term of linear stability in |PD93| . 
and recently in term of Lyapunov stability in jLMRj and [BCQlj by different 
ways. 

Theorem 4.5 The stability of the ring of N identical vortices depends on N 
and the latitude as follows: 

N=2 is Lyapunov stable at all latitudes; 

N=3 is Lyapunov stable at all latitudes; 

N=4 is Lyapunov stable if cos 2 #o > 1/3, and linearly unstable if the inequality 
is reversed; 

N=5 is Lyapunov stable if cos 2 8 'q > 1/2, and linearly unstable if the inequality 
is reversed; 
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Figure 5: Stability results for geostrophic Cn(R,p) relative equilibria. A relative 
equilibrium can be formally stable (S), linearly unstable (U), elliptic (E). The 
diagrams for N > 8 are similar to that for N = 11. 
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N=6 is Lyapunov stable if cos 2 9q > 4/5, and linearly unstable if the inequality 
is reversed; 

N>6 is always linearly unstable, 

where 9q is the co-latitude of the ring. 

The stability of Cn v (R,p) relative equilibria is given in the next theorem 
and illustrated by Figures El and [7| We assume that the polar vortex is at the 
North pole without loss of generality, and that the momentum of the relative 
equilibrium is non-zero. Let 0$ be the co-latitude of the ring. 

Theorem 4.6 (Laurent-Polz, Montaldi, Roberts (EMRJ) 

A C2v{R,P) relative equilibrium is Lyapunov stable if 

(1 + 2cos6 l o)[(l + cos6> ) 2 A + cos6» (2 + 3cos6> )] < 0. 
and linearly unstable if the inequality is reversed. 

A C^ v {R,p) relative equilibrium is Lyapunov stable if aX(X + n cos 9q)(X — Ai) < 
0, and spectrally unstable if and only if 

8aX > (nsin 2 9 + 4(n - 1) cos6l ) 2 . 

A Cnv{R,p) relative equilibrium with N > 4 

(i) is spectrally unstable if and only if 

X < X or 8aX > (N sin 2 8 a + 4(JV - l)cos6» ) 2 , 

(ii) is Lyapunov stable if 

A > A and aX(X + Ncos6 )(X - Ai) < 0, 

where 

a = (iVcos(9o-iV + 2)(l + cos6) ) 2 

Ai = (N- l)cos6> (^sin 2 6» +2(N- 1) cos <9 )/a 

A = (c A r-(iV-l)(l + cos 2 6lo))/(l + cose'o) 2 

J iV 2 /4 if N is even, 

CN ~ \(iV 2 -l)/4 if N is odd. 

Let x e be a Cmv{R) or Cjy v (R,p) relative equilibrium on a sphere rotating with 
angular velocity ft. By the theory of perturbations, there exists a neighborhood 
U of zero, U C M, such that for all f2 S U, x e has the stability of the corre- 
sponding relative equilibrium on the non-rotating sphere. Hence the stability 
results of the two previous theorem persist in a neighborhood of zero for f2. 
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Figure 6: Bifurcation diagrams for C2 V (R,p) and C3 V (R,p) relative equilibria. 
The circles represent the eigenvalues of the mode 1. A relative equilibrium can 
be Lyapunov stable (S), linearly unstable (U), elliptic (E). 




Figure 7: Bifurcation diagrams for Cn v (R, p) relative equilibria. The bifurcation 
diagrams for N > 7 are similar to that for N = 8, while those for N — 4, 6 are 
similar to that for N = 5. The circles represent the eigenvalues of the mode 1, 
while the cross represent those of the mode [N/2]. A relative equilibrium can 
be Lyapunov stable (S), linearly unstable (U), elliptic (E). 
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4.3 Conclusions 



From the results of that section, one can reasonably think that for a small 
angular velocity fl the relative equilibrium Cn v (R) with a small co-latitude 9q 
(that is the ring is close to the pole) is Lyapunov stable for N < 5, while it is 
linearly unstable for N > 8. The cases N — 6 and N — 7 are less clear since 
the stability for N = 6 depends on the parameter n in Proposition ^. 41 and the 
case N — 7 is stable or unstable depending on whether we are on the plane or 
the sphere. 

About the relative equilibrium Cnv(R,p) with a negative vorticity for the 
polar vortex and a small co-latitude 9q (that is the ring and the polar vortex 
vorticities of have opposite signs), one can reasonably think from the previous 
results that for a small angular velocity f2 the relative equilibrium is: 

• elliptic for N = 3, 

• elliptic for a weak polar vorticity and N = 4, 5, 6, 

• linearly unstable otherwise. 

5 The Southern Hemisphere Circulation 

It is well known from sailors that around Antarctica there are all year long 
strong winds coming from the West, the westerlies. The atmospheric circulation 
around Antarctica is indeed essentially zonal during all seasons of the year 
|RL53llS68) . the global circulation is a steady rotation of low pressure systems 
around the North-South axis. This circulation is commonly called The Southern 
Hemisphere Circulation by meteorologists. With the words of this paper, we can 
say that the atmospheric circulation is close to the relative equilibrium Cnv (R) 
formed of N low pressure systems. Since there exists a polar vortex (a high 
pressure system) at the South pole, the relevant relative equilibrium is merely 
Cnv(R,p)- The southern hemisphere circulation varies with seasons |L65I IL67j 
but also intraseasonally RL53, MFG91 : the westerlies are stronger in summer, 
and their latitude may fluctuate within a season, while the zonal character is 
always maintained. Polar high pressure systems, which have their origin in or 
near Antarctica, at times attain sufficient strength to break through the zonal 
flow to reach and reinforce a sub-tropical high pressure system (the anticyclonic 
system of Saint-Helen for example). 

Can point vortices explain the Southern Hemisphere Circulation? Indeed 
for sixty years some meteorologists or physicists study point vortices in order 
to understand motions of low/high pressure systems in the Earth's atmosphere 
[S43llMS71irF7Rlll)P98j . In [SHJ the linear stability of the geostrophic C 3 {R) 
is proved and according to the author, that result explains the presence of the 
three sub-tropical major high pressure systems on each hemisphere. For exam- 
ple in the southern hemisphere, we have the anticyclonic systems of Saint-Helen, 
of the Indian Ocean, and of the South Pacific, which persist. In the northern 
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hemisphere, the high pressure systems of the Acores is well known from euro- 
peans. In [F75 , a linear stability study of Cs v (R) in the /3-plane approximation 
is used to conclude that a ring of major high pressure systems is stable at sub- 
tropical latitude, while for low pressure systems it does not. However, a ring of 
low pressure systems is stable near the poles. The problem of those two papers 
is that sub-tropical latitudes are partially occupied by continents, and they do 
not take in account the continent-atmosphere interaction. 

In |DP98| , the authors relate their work on vortices to atmospheric blocking 
events. The advantage of the Southern Hemisphere Circulation is that it occurs 
around Antarctica, which is a continent with a disc shape, thus the SO(2) sym- 
metry is not broken. Moreover, the change continent-ocean on the Antarctica 
coast must influence the atmospheric circulation, but — again — this does not 
break the SO{2) symmetry. Note also that this region (40-90 South) is the only 
one large region of the Earth which has an SO(2) symmetry. 

From the previous section, one can conjecture that a ring of N low pressure 
systems is Lyapunov stable for N = 3, 4, 5; and a ring of N = 3, 4, 5, 6 low pres- 
sure systems together with a weak polar high pressure system is elliptic, while 
the whole system becomes unstable if the polar high gets stronger. This agrees 
with the observations of the beginning of this section. But first, we did not take 
into account the sub-tropical highs comparing the southern hemisphere circula- 
tion to Cn v (R,p) relative equilibria. Indeed the sub-tropical highs may balance 
the influence of the polar high, and therefore reinforce the stability of the ring 
of low pressure systems. Second, we did not take into account the external 
heating by the Sun which adds dissipation to the system, and we know that dis- 
sipation generically induces instability for elliptic relative equilibria. However, 
this instability may be compensated by other phenomena such as the presence 
of sub-tropical highs. With the external heating by the Sun, the vorticity is 
no longer a conserved quantity, and one prefers the potential vorticity which 
is almost conserved (see |DL93| ). The reader may have a look on the forecast 
for the southern atmospheric circulation on the web at the following address: 
http : //www. ecmwf . int/products/f orecasts/d/ charts/deterministic/world/. 
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Appendix A 



Let Bk be the sum 



sin(0 - cj)j. k ) 



~[^~ rcos((j) - 4>j t k) 

where cj)j.k = ffc + 2-7r(j — 1)/N, r = sin 6> sin 0fc/(l — cos 8 cos and is fixed. 
We have 

N k oo 
j — 1 m— 

for all \r\ < max{l/\ cos(<p — <pj,k)\,j = 1, • • • , Nk}- Since there exist ato, . ■ . , a m 
such that cos m (cj)-<f)j,k) = J2™=o a v cos(p(0 we have B fe = X)m=o 
where 

C m = ^2 a p ^ sin (^ ~ <t>j,k) cos(p(cf) - (f>j,k))- 
p=o i=i 

It follows from 53j=i sin(q((/) — (f)j,k)) = that C m = for all m, hence Bk = 
for all fc. 



Appendix B 

We give in this appendix the proof of the stability of the Cn(R,p) relative equi- 
libria on a non-rotating plane, that is the proof of the following theorem: 
Let En = +1 if N is even, and = if N is odd. 

The Cn(R,p) relative equilibria with N > 4 on a non-rotating plane are Lya- 
punov stable if max(0, (iV 2 - 8iV + 7 + £jv)/16) < A < (N - l) 2 /4, elliptic if 
(N 2 — 8N + 7 + £jv)/16 < A < 0, and linearly unstable if the previous condi- 
tions are not satisfied. The Cz(R,p) relative equilibria are Lyapunov stable if 
< A < 1, elliptic if X < 0, and linearly unstable if X > 1. 

The proof is quite long though we omitted lots of details. We first recall some 
well-known facts about the dynamics of planar point vortices. 

Let x e be a Cn(R,p) relative equilibrium and A be the vorticity of the central 
vortex. The equations of motion N + 1 planar point vortices are |Ar82j : 

_ 1 ^ A fc 



3 2m £-i . Zj - z k 

where Zj is a complex number representing the position of the j-vortex (we iden- 
tified the plane with C), and Xj is the vorticity of the j- vortex. The Hamiltonian 
for this system is 



H = - — > Xi A,Tn | z t - *j | 
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and the symmetry group of the vector field Xh is SE(2) (which is not compact). 
In the particular case of N identical vortices together with an additional vortex, 
the Hamiltonian is E(2) x S^r-invariant. Let Zj = pj exp(i<f>j), j = 1, . . . , N and 
zjv + i = x + iy, the dynamical system in (fij, (pj, x, y) variables is: 

w> = iM ' * = 

After identifying SE(2) with C XI 5 1 and so se(2)* with Cxi, the momentum 
map of the system is 

N+l JV+l 

i=i i=i 

In a frame such that z^ + i = 0, a Cpj(R,p) configuration has coordinates pj — 
R,(pj — 2iTj/N,j = l...N,x = 0,y = 0, and is a relative equilibrium with 
angular velocity £ = (N — 1 + 2A)/4i? 2 (where H was normalized by Air). It 
is straightforward to see that R 2 H^ does not depend on R, hence so does the 
stability, and we can set R = 1 . 

The recent results of Patrick, Roberts and Wulff PRW02 generalize the 
Energy-Momentum method to Hamiltonian systems with a non-compact sym- 
metry group. However, we will apply here the classical Energy-Momentum 
method with the compact subgoup 50(2), and thus forget the translational 
symmetries. The momentum map coming from the rotational symmetries is 
J = X(x 2 + y 2 )/2 + ~^2j Xjp 2 /2. Since the coadgoint action of 50(2) is trivial, 
one has 50(2) M = 50(2). The symplectic slice at a Cn(R,p) relative equilib- 
rium is therefore 



Af = Span 



( N N 



Let G = 0(2) x Sn- The linear map associated to d 2 H^\j^(x e ) (see Intro- 
duction) is equivariant under both symplectic and anti-symplectic elements of 
G Xe , while L/v is equivariant under only the symplectic ones (see |LP02| and 
LMR for a detailed account on the sphere). We then perform an isotypic 
decomposition and find that the symmetry adapted basis of M is 

(«f 5 ,$\Ty^\4\T X ,{ af , af , fl> , ff$> | 2 < / < [n/2] }) 

where 

af+iffl = J2s=i e M^£s/N)5p s 
<*f+i$ = Ef=i exp(2in£ S /N)S ( p s 

5x = --^=5x 

Sy = ~^=5y. 



20 



In this basis, the matrix d 2 H^\j^{x e ) (resp. Ljv) block diagonalizes in 2 x 2 
blocks with two 3x3 blocks (resp. 4x4 blocks with one 6x6 block, plus a 
2x2 block if N is even) . 

We first study the Lyapunov stability. Some simple calculations show that 



d 2 H i (x e ) 




= d 2 H ( (x e 


d 2 H i (x e ) 




= d 2 H i (x e 


^H^Xe) 




= d 2 H i {x e 


d 2 H i (x e ) 




= d 2 H i (x e 
= d 2 H e {x e 


d 2 H i (x e ) 




= d 2 H i {x e 



a g ,a e ) 



«e >°Q, ) = 



atg ) jSx) = NXy/N/2 
8x,5x) = ATfA/2. 



Hence it follows from the block diagonalization that 

d 2 H i \ M {x e ) = diag{A,A,D) 
where D = diag({\^ , Xf , Xf\ \2<£< [n/2]}), 



.4 





a N 







N£X/2 



and 



A 



(<') 



Oat 



d 2 H i (x e 
d 2 H^x e 



= NXJN/2. 



Note that £> exists only if AT > 4. 
Thanks to the formula H75 



JV-l 

E 



cos(2n£j/N) 
sin 2 (7rj/7V) 



1) - 2£(N-£) 



wc find after some lengthy computations that = N£{N — l)/2 and 



The eigenvalues A^f" 1 are all positive, thus D is definite if — (£— 1)(N — £ — 1) + 
N - 1 + 4A > for alU = 2, . . . , [AT/2], that is if A > (([AT/2] - 1)(AT - [AT/2] - 
1) - N + l)/4 which corresponds to A > (N 2 - + 8)/16 for AT even and to 
X> (N 2 ~8N + 7) /16 for AT odd. 

The relative equilibrium is therefore Lyapunov stable if A is positive definite, 
that is if the three following subdeterminants are positive: 



^[-(£-l)(N-£- 



1)+N - 1 + 4A]. 



A 



(i) 



detA 
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Since A 



1 1 --N(N- l)/2, A^ 
det,4 = -TV 3 A 



TV(4A + N 
A-^ 



l)/2 and 



it follows that A is positive definite if and only if < A < (TV — l) 2 /4. 

We proved therefore that Cn{R,p) relative equilibria with TV > 4 are Lya- 
punov stable modulo 50(2) if max(0, {N 2 -8N + 7 + e N )/16) < A< (TV-l) 2 /4 
where En — (resp. En — 1) for TV odd (resp. even). 

We then study the linear stability of the relative equilibrium. It follows 
from the block diagonalization of L/v that L/v — diag {Al, {At | 2 < £ < [n/2]}) 
where 



At, = 



I 











m/2 


a N /X 


a N /X 













-a N 


















-ajv 


iM 









-aw/A 


~a N /X 













a N 





A e 











V 




At) 
x 















and 



/ 





A 
















\W 
A e 


>(<) 















for 



while it is 







= N/2 (TVeven) 



(the blocks are given up to a factor). 

After some calculations, we obtain that the eigenvalues of are 



4 & , ± 



(N-iy 



x 



±WAW|2<K[n/2] 



We have therefore some double eigenvalues, but the Jordan forms of the blocks 
are semi-simple. It follows from the discussion on Lyapunov stability that the 

eigenvalues ±i 



X^p A^ are all purely imaginary if and only if A > (TV 2 — 8TV + 
7 + £jv)/16. The Theorem follows for TV > 4. 

For TV — 3, the Lyapunov (resp. linear) stability is governed only by A (resp. 
Al). Hence Cs(R,p) is Lyapunov stable ifO < A < (TV — l) 2 /4 and linearly 
unstable if and only if A > (TV - l) 2 /4. 
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